Controlling transistor threshold voltages using molecular dipoles 
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We develop a theoretical model for how organic molecules can control the electronic and transport 
properties of an underlying transistor channel to whose surface they are chemically bonded. The 
influence arises from a combination of long-ranged dipolar electrostatics due to the molecular head- 
groups, as well as short-ranged charge transfer and interfacial dipole driven by equilibrium band- 
aUgnment between the molecular backbone and the reconstructed semiconductor surface atoms. 



Inorganic semiconductors have traditionally domi- 
nated as the material players in the electronics indus- 
try. While their organic counterparts have been studied 
extensively as alternate channel materials [l[, the devel- 
opment of a stand-alone molecular electronics technology 
has been stymied by the inordinate difficulty of contact- 
ing small molecules reproducibly, their insufficient mobil- 
ities, large RC constants and poor gateability. Perhaps 
a more promising approach is to envisage hybrid organo- 
semiconductor devices, combining the established infras- 
tructure of the semiconductor integrated industry with 
the 'bottom-up' self-assembly and chemical tunability of 
molecular monolayers. A particularly interesting possi- 
bility is to use organic molecules to control the surface 
properties of deeply scaled, backgated silicon transistors, 
by tying up deleterious surface states and charge transfer 
'doping'. In addition, monitoring the transistor dynam- 
ics, such as the shift in its threshold voltage, can be used 
to detect a single molecule adsorption event [2| . It is thus 
critical to properly understand the physical factors that 
determine how a molecule controls a transistor, and how 
the transistor, in turn, senses the molecule. 

In this paper, we develop a quantitative theory for the 
threshold voltage control of low-doped silicon channels by 
surface bonded organic monolayers with varying dipole 
moments. We focus in particular on a recent series of 
experiments [4| that involved the grafting of molecular 
monolayers atop oxide free H-passivated silicon surfaces. 
The choice of the molecules followed an important logic. 
An identical set of molecules was used, with the excep- 
tion of one substituent group. This allowed a systematic 
study of the effect of the molecules on the electrical prop- 
erties of the device. While organic molecules attached 
to semiconductor surfaces have been studied extensively 
[51, lo, lZI, albeit phenomenologically, our principal chal- 
lenge is to develop a quantitative, 'first principles' model 
that combines atomistic charge-transfer processes with 
macroscale electrostatics. We employ Density Functional 
Theory (DFT) to extract the molecular adsorption geom- 
etry, interfacial dipole and band- alignment at the atom- 



istically reconstructed silicon surface. These quantities 
are then incorporated as inputs into a macroscopic Pois- 
son solver to compute the band-bending in the transis- 
tor channel. The calculated threshold voltage shifts and 
band-alignments are in excellent agreement with experi- 
ments for multiple molecules (Table I) [J]. 

Given the low doping and likely statistical fluctuations 
in data, as well as uncertainties in the experiment (for 
instance, in the number of bonded molecules), there is 
undoubtedly some wiggle room for theory that may make 
this correspondence fortuitous. However, the main trend 
in the experiment seems quite robust, namely, a clear 
monotonic dependence of the threshold voltage on the 
dipole moment of the headgroup. This is borne out quite 
well by our model. In addition, we demonstrate the accu- 
racy of our model by comparing our computed interme- 
diate ingredients with complementary experimental data 
obtained on the same set of samples [8| . 
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FIG. 1: Schematic side-view representation (not to scale) of 
the pseudo-MOSFET under investigation [3]. The molecules 
were grafted between the source and drain electrodes. Vs, 
Vd, and Vg refer to the bias applied on the source, drain, and 
gate, respectively. 

Molecule free device. The devices we are studying were 
fabricated using a silicon-on insulator (SOI) wafer [J]. 
A back gated structure was used to allow easy grafting 
of a molecular monolayer on the silicon device top layer 



(FigH]). The molecules were grafted on H-passivated sili- 
con surfaces using diazonium based chemistry (molecules 
lose the diazoniumgroup to form radicals before attach- 
ing to the surfaceQ). In the fabricated MOSFET, the 
bulk p-Si substrate (handle) was biased at Vg using the 
gold back contact and acts as a gate terminal. This in 
turn induces a conduction channel at the upper interface 
of the buried oxide, used as a gate dielectric layer. The 
conducting channel corresponds to the formation of an 
accumulation layer (p-channel). In practice, the drain 
current voltage characteristic presents a linear behav- 
ior at low drain voltages Vd, and tends to saturate as 
the drain voltage approaches the voltage difference be- 
tween the gate bias and the threshold voltage. In the 
diffusive transport regime, the current- volt age character- 
istics are satisfactorily described by the square law theory 
0(Figll 
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where W and L are the channel width and length, Co 
is the gate capacitance and /i is the channel mobility. In 
this case, an accumulation channel is activated when both 
the gate and the drain are negatively biased. Because 
the transistor body is the nearly intrinsic p-Si layer, the 
channel is assumed to be completely accumulated. 
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FIG. 2: (a) Experimentally obtained I-V characteristics [3] of 
the fabricated MOSFET devices. Inset (a) displays the typ- 
ical transfer characteristics of the devices under test. Inset 
(/?) shows the typical output characteristics of the devices 
under test. (b) Theoretically calculated I-V characteristics of 
the MOSFET device. 

An accumulation mode transistor is not expected to 
have a prominent turn-on as one expects for an inver- 
sion channel transistor. The onset voltage is associ- 
ated with pulling the contact Fermi energies below the 
valence band-edge near the oxide interface. The large 
threshold voltage values of the bare MOSFET (- -5.6 V) 
could be indicative of traps at the oxide-channel interface. 
This is further confirmed by the presence of prominent 
hysteretic signatures that suggest charge trapping pro- 
cesses. We used two approaches to estimate the areal 
density of inter facial traps in the control devices. The 
first assumes that the threshold voltage of the control 
devices is due to trapped charge on one side of a capac- 
itor, with a capacitance/area equal to Cqx- Using this 



method, the density of interfacial traps was calculated to 
be Ni = Qp/q = \AVt\ x Cc/q = 1.2 x lO^Vcm^ The 
second approach estimates the same quantity using the 
observed subthreshold swing S [uj . The traps would de- 
grade the subthreshold current by taking away some of 
the mobile charges that would otherwise have contributed 
to current flow. The interfacial trap density is calculated 
to be Ni = SCG\og^oe/{kT/q) = 1.0 x lO^Vcm^. The 
consistency of the two estimates strengthens our argu- 
ments as to the origin of the threshold voltage and sub- 
threshold slopes observed. 

Molecular attachment. Quantum-chemical calcula- 
tions were carried out using DFT as implemented in the 
Vienna Ab-initio Simulation Package (VASP) [H, [l3|. 
Plane wave basis sets were used to represent the atoms 
and the PW91 form of the GGA functional was used to 
carry out gradient corrected calculations. The optimized 
structures for the isolated gas-phase molecules under con- 
sideration, dimethylaminobenzene, aniline and nitroben- 
zene and the corresponding radical forms are shown 
in FiglH These structures were obtained using spin- 
polarized calculations, and are in agreement with pub- 
lished experiments and theoretical computations [ij, [l5| . 

The Si(lOO) slab was modeled using 9 atomic layers 
and hydrogen atoms were attached on the lowest layer of 
the slab to passivate the dangling bonds on the lower slab 
surface. Geometry optimizations were carried out while 
constraining the lowest 4 layers of the Si(lOO) slab at their 
bulk positions including the hydrogen atoms attached on 
the lowest atomic layer. We started with a dimerized 
Si(lOO) structure as our initial guess (with broken sym- 
metry), and recovered the expected 2x1 reconstructed 
structure. The dipole moments and vibrational modes 
were calculated within VASP using single point calcula- 
tions. 

Once our isolated geometries were thoroughly bench- 
marked, we moved on to the adsorption geometry. The 
radical forms of the molecules we studied, specifically, 
dimethylaminobenzene, aniline and nitrobenzene, were 
adsorbed on the asymmetrically dimerized Si(lOO) slab 
via a C-Si covalent bond. We considered two possible 
adsorption sites on the Si(lOO) slab, named Tl and T2. 
In all the cases considered for atop adsorption of radicals 
on the surface, we find that the radicals stand normal 
to the surface. A positive dipole, with the positive pole 
directed toward the semiconductor surface, is induced by 
nitrobenzene, while a negative dipole is induced by ani- 
line and dimethylaminobenzene. Table I shows the com- 
puted dipole moments. 

Threshold Voltage Shift. The presence of the molecu- 
lar layer modifies the silicon surface work function, which 
in turn affects the device properties. The work function 
(WF) is the minimum energy required for an electron to 
escape into vacuum from the Fermi level {Ep) of the ma- 
terial, and is determined by (i) the electron affinity (EA), 
the energy required to excite an electron from the bottom 




FIG. 3: Optimized geometries for the (a) gas-phase molecule 
aniline(arrow points in the direction of dipole moment, from 
positive to negative charge). (b) Optimized structure for pas- 
sivated Si(lOO) slab with symmetrically dimerized surface. (c) 
Optimized structure for aniline radical adsorbed on passivated 
Si(100)-(2xl) slab atop site Tl( top atom of the dimer). 



of the conduction band (CB) at the surface to the local 
vacuum level; (ii) the band bending (BB), the electrical 
potential difference between the surface and the electri- 
cally neutral semiconductor bulk; and (iii) the energy dif- 
ference between the bulk CB and the Fermi level. AWF 
can be due to AEA, ABB, or both. We computed the 
total shift directly from the shift in the estimated highest 
occupied molecular orbital (HOMO) level of the molecule 
before and after adsorption. This was accomplished by 
comparing the projected density of states on the sili- 
con atoms, which is maximized for the HOMO levels of 
benzene-based molecules [16]. The computed shifts agree 
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FIG. 4: Comparison of calculated values and experimental 
UPS/IPES/XPS data. Experimental data are adapted from 



very well with experimental data, as seen from the last 
two columns of Table I, as well as comparisons with ultra- 
violet photoelectron spectroscopy (UPS) and inverse pho- 
toemission spectroscopy (IPES) data combined with X- 
ray photoelectron spectroscopy (XPS) data (FigH]) The 
shift shows a monotonic dependence on the dipole mo- 



ment of the headgroup, which is quite encouraging. Sur- 
prisingly, however, there is an additional shift relative to 
the H-passivated control, which is unexpected (thus, the 
threshold voltage shifts of the nitro and amino compo- 
nents do not straddle that of the control, but lie on the 
same side of it). We will now attempt to deconstruct 
the computed shifts to understand their physical origins. 
Dipolar Contribution. The electron affinity of a given sur- 
face is directly affected by a surface dipole, which creates 
an electrical potential drop across the grafted molecu- 
lar film depending on the dipole moment. Simply put, 
a positive headgroup such as NH2 acts as a gate that 
pulls negative charges to the surface, lowering its local 
electronic levels, while NO2 pushes them away (Fig|5]). 
The resulting band-bending at the surface propagates to- 
wards the opposite end where the oxide sits, penetrating 
up to a Debye length L^ = Y^e^eo/c^T/A^^g^ , where e^ 
is the dielectric constant of the channel, eo is the permit- 
tivity of free space, T is the temperature and Na is the 
acceptor dopant density. 
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FIG. 5: Adsorbed layer of (a) Aniline (Blue line) and 
(b) Nitrobenzene (Brown line) molecules, creating dipole mo- 
ments on the surface, (c) Energy Band Diagram for the sys- 
tem being studied. (d) Pictorial representation of change in 
band bending due to charge transfer between molecular mono- 
layer and substrate. 

The purely dipolar headgroup gives a shift in potential 
given by 

¥dip = {N^oil^dip cos 6>/e^)e-^-/^^ (2) 

where Nmoi gives the areal surface density of the attached 



TABLE I: Data of computed dipoles of the molecule-semiconductor system, shifts due to dipole {(j)dip) and charge transfer 
((/)bb), computed and experimental changes [4] in the threshold voltage of the devices. The control consists of an H-passivated 
Si(lOO) surface. 



Systems fidip ^(l^dip (^066 AVt(V^) AVt(V^) 
(dB) (V) (V) Theory Expt. 

Control 0.03 0.0058 0.0058 

Nitro 2.17 0.22 -0.725 -0.525 -0.55 

Aniline -1.44 -0.14 -0.96 -1.1 -1.25 

Dimethylaminobenzene -1.65 -0.28 -1.12 -1.4 -1.55 



molecules, fidip is the static dipole moment of each molec- 
ular head group oriented at an average angle relative to 
the normal to the surface, e^ is the molecular dielectric 
constant, and tsi is the thickness of the channel silicon. 
This contribution is opposite for opposite dipolar signs. 
Significantly, dipolar gating at the top surface nontriv- 
ially influences the threshold voltage at the backgate ow- 
ing to the large Debye length associated with low doping 
in the channel. For a doping level of IO^^ctti"^ the De- 
bye length Ld ~ 1.25/im, so that the exponential trans- 
fer factor is around 0.7 (implying that seventy per cent 
of band-bendings at the surface are transmitted to the 
bottom end). 

Charge- Transfer Barrier. In addition to the dipo- 
lar contribution, the molecular monolayers also trans- 
fer spectral weight to the semiconductor surface through 
bonding. Because the mobile charges are primarily near 
the bottom of the channel, this does not add signifi- 
cant additional resistance, but it does influence the local 
band-bending through the charge-transfer dipole at the 
molecule-semiconductor interface. Based on our com- 
puted shifts and the dipolar contributions separately 
extracted, we can compute this component for each 
molecule. To zeroth order, this represents the charge 
transfer due to the workfunction difference between the 
headgroup-free molecular backbone (the bare benzene 
ring) and the silicon surface. This means that the as- 
sociated band-bending is in the same direction for all 
the headgroups (FiglSb), and relatively weakly varying 
with dipole moment (column 4 of Table I). There is, 
however, some variation, indicating that the separation 
into a purely dipolar part and a head-group indepen- 
dent backbone part is not strictly feasible, owing to the 
short length of the benzene rings and the resulting hy- 
bridization between the ring states and the headgroups. 
The electron donating ring electrons push charge to- 
wards the p-Si substrates and become partly positively 
charged, leading to an increase in the positive surface 
charge density and a lowering of the local silicon sur- 
face levels (FiglSb). The donating capacity, however, 
is enhanced (diminished) by the electronegativity (pos- 
it ivity) of the headgroups. Methylamine has the high- 
est electron-donation capability, and hence provides the 
highest change in charge-transfer induced band bending 
^^66- While we do not independently compute ^^55, one 



can roughly estimate this by calculating the charge neu- 
trality level EcNL of the headgroup molecular backbone 
(CNL is defined as the energy to which electrons need 
to fill the molecule to keep it electrically neutral - fre- 
quently this is replaced by the LUMO level H). The 
work-function shift of the molecule is then given by 



{Ep — Ecnl)^' 



-tch/Lo 



^ 1 + l/UoDo ^ ^ 

where Uq is the single-electron charging energy of the 
molecule and I^o is its density of states near the Fermi 
energy. The barrier height Ep — Ecnl is of the order 
of 1 V for aromatic molecules on silicon [7[, giving us 
the correct observed order of magnitude. The term in 
the denominator includes the Coulomb cost for charg- 
ing up the molecule. In this contribution, we have stud- 
ied the effects, both qualitatively and quantitatively, of 
grafting molecular monolayers on the top surface of back- 
gated MOSFETS. Such a scheme of using molecules on 
surfaces to tailor device properties is becoming increas- 
ingly attractive as devices get smaller, and surface-to- 
volume ratios increase vastly. The study of surface mod- 
ifications could also have practical implications for sen- 
sors, specifically chemically sensitive field-effect transis- 
tors (CHEMFETs)[i3. A radicallydifferent principle of 
operation that we are working on [l8| involves a stronger 
bonding between the molecule and the channel, involving 
actual wavefunction overlap, leading to quantum scat- 
tering by the molecular traps that creates characteristic 
fingerprints when scanned with a back gate. 
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